Quantum stochastic models of two-level atoms and electromagnetic cross sections. 
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Quantum stochastic differential equations have been used to describe the dynamics of an atom 
interacting with the electromagnetic field via absorption/emission processes. Here, by using the full 
quantum stochastic Schrodinger equation proposed by Hudson and Parthasarathy fifteen years ago, 
we show that such models can be generalized to include other processes into the interaction. In the 
case of a two-level atom we construct a model in which the interaction with the field is due either to 
absorption/emission processes either to direct scattering processes, which simulate the interaction 
due to virtual transitions to the levels which have been eliminated from the description. 

To see the effects of the new terms, we study various types of cross sections for the scattering 
of monochromatic coherent light. We obtain formulas giving the total, the elastic and the inelas- 
tic cross sections as functions of the frequency and the intensity of the stimulating laser and the 
fluorescence spectrum as a function also of the frequency of the scattered light. The total cross 
section, as a function of the frequency of the stimulating laser, can present not only a Lorentzian 
shape, but the full variety of Fano profiles; intensity dependent widths and shifts are obtained. 
The fluorescence spectrum can present complicated shapes, according to the values of the various 
parameters; when the direct scattering is not important the usual symmetric triplet structure of the 
MoUow spectrum appears (for high intensity of the stimulating laser), while a strong contribution 
of the direct scattering process can distort such a triplet structure or can even make it disappear. 
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I. INTRODUCTION 

Quantum stochastic calculus (QSC) [|l]-D, a noncom- 
mutative analog of the classical Ito's stochastic calculus, 
revealed to be a powerful tool to construct mathematical 
models of quantum optical systems and to de- 

velop a theory of photon detection |13|-|16[. Just at the 
beginning of QSC, Hudson and Parthasarathy proposed 
a quantum stochastic Schrodinger equation for quantum 
open systems [QU]. Such an equation has found appli- 
cations in quantum optics, but not in its full generality 
|6|,P^|J^. It has been used to give, at least approximately, 
the dynamics of photoemissive sources such as an atom 
absorbing and emitting light, or matter in an optical 
cavity, which exchanges light with the surrounding free 
space. But in these cases the possibility of introducing 
the so called gauge (or number) process in the dynam- 
ical equation has not been considered; roughly speak- 
ing, the gauge process is a quadratic expression in the 
field operators which preserves the number of quanta, 
but changes their wave functions. In this paper we want 
to show, in the case of the simplest photoemissive source, 
namely a two-level atom stimulated by a laser, how the 
full Hudson-Parthasarathy equation allows to describe in 
a consistent way the scattering of the light by the atom 
not only through the absorption/emission channel, but 
also through another process which can be called "direct 



scattering" , which can be included in the interaction via a 
term containing the gauge process. When the atom is ap- 
proximated by a two-level system, the introduction of an 
interaction term which preserves the number of photons 
allows to simulate also the scattering processes involving 
virtual transitions to states different from the two ones 
responsible of the real absorption/emission process. 



So, a first aim is to show how the full Hudson- 
Parthasarathy equation is able to give a reasonable and 
rich model for the dynamics of an atom interacting with 
the electromagnetic field. A second one will be the study 
of the elastic, inelastic and total cross sections for the 
scattering of monochromatic coherent light by the atom. 
The resulting line-shapes are very interesting. For in- 
stance, the dependence of the total cross section on the 
frequency of the stimulating laser can present not only 
a Lorentzian shape, but the full variety of Fano profiles 
(]l7t, Q pp. 61-63). Moreover, the dependence of the 
line shape on the intensity of the stimulating laser is 
computed and power broadening and intensity dependent 
shifts are found. The study of the inelastic cross section, 
instead, shows possible modifications to the known triplet 
structure of the fluorescence spectrum [ p^ . Some prelim- 
inary results on the total cross sections where reported 
in 
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Fock space and QSC 

Let us recall some notions of QSC and the Hudson- 
Parthasarathy equation; this is just to fix our notations, 
while for the proper mathematical definitions and the 
rules of QSC we refer to the book by Parthasarathy 
We denote by ^ = !F{X) the Boson Fock space over 
the "one-particle space" X = Z (» L'^{R+) L'^{R+; Z), 
where Z is another separable complex Hilbert space. A 
vector / in A" is a function from R_|_ into Z; we fix a 
c.o.n.s. {ci, i > 1} in Z and we denote by 



Mt) = {ej\f{t)) 



(1.1) 



the components of a vector f{t) in Z. The Fock space T 
is spanned by the exponential vectors E(f), whose com- 
ponents in the 0, 1, . . . , A:, . . . particle spaces are 

E{f) - (l, /, (2!)-i/2/ ^ /, . . . , . . .) , (1.2) 

f €z X] the inner product between two exponential vec- 
tors is given by 



{E{g)\E{f)) = e^p (ffl/) = exp 

r + co 



{9m{t))dt 



exp 



(1.3) 



where an overline means complex conjugation, and we 
get normalized vectors by defining 



K/)=exp(--i||/|p)i?(/). 



(1.4) 



The annihilation, creation and gauge (or number) pro- 
cesses are defined by 



Ajit)E{f)^ / f,is)dsEif), 
Jo 

(Eig)\Alit)Eif)) = f^ds{E{g)\E{f)), (1.5) 
Jo 

{E{g)\k,,{t)E{f)) = f ^f,{.s)d.s{Eig)\E{f)). 



Eqs. (O) allow to write formally 



aj{s) ds . 



Alit)^ / a]{s)ds, 
Jo 

Kij{t) = / a\{s)aj{s)ds , 



(1.6) 



where aj(t), aj(i) are usual Bose fields, satisfying the 
canonical commutation rules 

K-(t),a,(s)] =0, [a,[t),a\[s)]^5,,5[t~s), (1.7) 



and whose coherent vectors are the normalized exponen- 
tial vectors: 



a,{t)eU)^f,{t)e{f). 



(1.8) 



In particular the vector e(0) = E(Q) is the Fock vacuum. 

The Bose fields introduced here represent a good 
approximation of the electromagnetic field in the 
so called quasi-monochromatic paraxial approximation 
pl| , p^ . Now, T is interpreted as the Hilbert space of 
the electromagnetic field; ^j(i) creates a photon with 
state Cj in the time interval [0,t], Aj{t) annihilates it, 
^jjit) is the selfadjoint operator representing the num- 
ber of photons with state Cj in the time interval [0, t\ 
and 



7V(t)=EA,,(i) 



(1.9) 



is the observable "total number of photons entering the 
system up to time . Moreover, in the approximation we 
are considering, the fields behave as monodimensional 
waves, so that a change of position is equivalent to a 
change of time and viceversa. If we forget polarization, 
the one-particle space Z has to contain only the degrees 
of freedom linked to the direction of propagation 1 22 1 , so 
that we can take 



Z = L2(T, sinede'd^) , 

T = {0 < 6* < TT, < < 27r} ; 



(1.10) 



the angular coordinates (0, </>) represent the direction of 
propagation. Now, a vector / in the one-particle space 
X can be identified with a function f{d^(j)^t) such that 

/+°°dt/;d0 sin0/o'"d0im0,t)i' < -i-oo. 

In QSC integrals of "Ito type" with respect to dAj{t), 
dAj(t), dh-ijit) are defined. The main practical rules to 
manipulate "Ito differentials" are the facts that dAjit), 
dA^jit), dkij (t) commute with anything contain the fields 
only up to time t and that the products of the fundamen- 
tal differentials satisfy 

dAj{t)dAl{t)^Sj,dt, 
dA,{t)dAk^{t) = 5,kdA,{t), 

dAj,{t)dAl{t)^6,kdAl{t), (1.11) 
dAj,(t) dAik{t)^S,i dAjfe(i), 

dAlit) dAj{t) = dAfe,(i) dA,{t) = dAlit) dA,,(t) = ; 

all the products of dAj (t) , dAj (t) or dA^ j (t) with dt van- 
ish. 
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The evolution equation 



Let 7i be a separable complex Hilbert space (the 



system space) and let Ri 



> 1, Sr 



> 



1, 

H be bounded operators in Tl such that = H, 
J2i R\Ri is strongly convergent to a bounded operator, 
and ^ - ■ Sij®\ei) {ej \ = S e U{T-L®Z) (unitary operators 
inH.® Z); we set also 



K = H 



(1.12) 



Then ( ||] Theor. 27.8 p. 228) there exists a unique 
unitary operator- valued adapted process U{t) satisfying 
L/(0) = 1 and 



dL/(t) = |^i?,d4(t)+^(5. 

3 ij 



d^j W ~^Kdt\ U{t) . (1.13) 

The operator Ut will be the evolution operator for the 
atom-field system, in the interaction picture with respect 
to the free dynamics of the field. In order to describe a 
two- level atom, we take H — C'^; then, to fix the model, 
we have to determine the atomic operators iJ, Ri, Sij 
on the basis of physical considerations. In the next sec- 
tion we shall require: (a) the existence of a ground state 
to which the atom decays by emitting at most one pho- 
ton when it is not stimulated, (b) a balance equation 
[Eq. ( 2.16| )] between the numbers of ingoing and outgo- 
ing photons when there is some coherent source. This 
suffices to determine the structure of the atomic opera- 
tors [Eqs. (p^, <^J%]. 



Contents 

As said before, in Section |l| we fix the model by phys- 
ical considerations; here, a central role is played by a 
balance equation saying that the mean number of outgo- 
ing photons plus the mean number of photons stored in 
the atom is equal to the mean number of ingoing pho- 
tons. In Section lU we consider the case of a spherically 



light, we obtain the master equation which gives the time 
evolution of the reduced atomic density matrix and we 
study the large-time behavior of its solutions. In Section 



IV 



we study the differential (with respect to the angle) 
and total cross sections for the scattering of laser light by 
the atom, as a function of the frequency of the stimulat- 
ing laser; in this section such cross sections are obtained 
from the direct detection scheme. In Section 0, starting 
from the balanced heterodyne detection scheme, we ob- 
tain the power spectrum of the fluorescence light and the 
elastic and inelastic cross sections. Section VI is devoted 



to a discussion of the main features of the integral cross 
sections and of the power spectrum. 



II. THE MODEL AND THE BALANCE 
EQUATION FOR THE NUMBER OF PHOTONS 

First of all we want a model for an atom stimulated by 
a laser (coherent light, not necessarily monochromatic); 
this means to choose as initial state 5'(^, f) € TL ® a. 
generic state for the atom and a coherent vector for the 
field im, i.e. 



*(e,/) = e®e(/), 

^&H, ||e|| = l, feL^{ 



(2.1) 



symmetric atom stimulated by monochromatic coherent 



Then, the atomic reduced statistical operator p{t; ^, /) 
is defined by the partial trace over the Fock space 

p(i;e,/)=Tr,{C/(i)|vI'(e,/))(vI'(^,/)|{7(t)t}. (2.2) 

Moreover, the quantity 

{N{t))f = {Um{i, f)\N{t) Um{^, /)) (2.3) 

represents the mean number of photons up to time t, 
after the interaction with the atom, while 

{N{t))) = {^{tf)\N{t)-^{iJ))= f\\f{s)fds (2.4) 

Jo 

is the same quantity before such an interaction |l^]; we 
can also say that Eq. ( p.4| ) gives the mean number of ingo- 
ing photons entering the system in the time interval [0, t] 
and that Eq. (2.3) gives the mean number of outgoing 
photons leaving the system in the same time interval. 



Proposition 1 The reduced statistical operator p(t; ^, /) satisfies the master equation 

d 



dt 



p{t-tf)^C{f{t))[p{t-i,f)], 



where 



L{f{t))[p\ = -i [H{f{t)) , p] + ([i?^.(/(t))p, R,{f{t))^] + [R,{f{t)) , pR,{f{t))^ 



(2.5) 



(2.6a) 
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i 



(2.6b) 



H{f{t)) ^H- '-J2(R]S,J.,{t) - Mt) S/,R,) ; 



(2.6c) 



moreover, we 



have 



{N{t))f = Tr„{^i?,(/(s))^i?,(/(s))p(s;^,/)}ds 



(2.7) 



Proof. By using the rules of QSC, it is possible to differ- 
entiate {N{t))f and {U{t)-^{i,f)\aU{t)-^{iJ)), where 
a is a generic system operator. Then, one gets the re- 
sults by recalling that the increments of the field op- 
erators commute with U{t) and that dAj{t)'^{^, f) = 
fj{t)dt'ii{^,f) and by using the definition of p{t;^,f) 
given in Eq. (p^). □ 

In order to formulate physical requirements, let us start 
by considering the case when no photon is injected into 
the system, i.e. / = 0. In these conditions it is natu- 
ral to ask that the atom can emit at most one photon; 
moreover, we require the existence of a unique equilib- 
rium state which we denote by pg (it will be the ground 
state). We take as canonical basis |— )} in Ti an 

orthonormal basis which diagonalises pg, so that we can 
write 



Pg = P- 



Pg = PP+ + {I - P)P- 



(2.8) 



for some p in [0, 1], where P± are the orthogonal projec- 
tions over the vectors |±). We shall use also the Pauli 
matrices 



1 

-1 

1 




-i 

1 



1 



(2.9a) 
(2.9b) 



by which the two orthogonal projections P± can be writ- 
ten as 

P+ = ^(1 + <yz) = <y+<y- , P- = ^(i - ff.) = <y-<y+ ■ 

(2.9c) 



Proposition 2 We require 

{N{t)) /=o < 1 , Ve, Vt, 



P{t] 0) Pg . 



(2.10a) 
(2.10b) 



Then, apart from an exchange of roles between the two 
states \+) and I—), we obtain 



a e Z , a 7^ 0. 



(2.11) 



(2.12a) 



(2.12b) 



Vice versa , Eqs. \2.12(\ ) and ( \2.12l\ ) imply Eqs. ( ^.1(\ ) 
and i2.1l\). 



Proof. By Eqs. (^ and ( p.lOaj) , (iV(i))/ is 
a bounded and non decreasing fun ction of t , so 
limt^+oo{N{t))f exists; then, Eqs. ( ^ and ( |2.10b| ) give 

J2j Tr„ ^RjRjPg^ — 0. By the cychc property of the 

trace and the positivity of pg and of RjPgRj, we get that 
this condition is equivalent to Rj pg = 0, Vj. 

Now, let us set Rj — Xji + HjCTz + + ctjO"- (every 
operator on can be written in this way). Then, Eq. 
(O) and RjPg = give p{xj+yj) = 0, {l-p){xj-yj) = 
0, (1 — p)zj = 0, paj ~ 0. For p g (0, 1) this system of 
equations gives Rj — 0, which is not acceptable because 
in this case the equilibrium state is not unique. For p = 
we get Xj = yj and Zj =0; we need also 



El 



^0 



(2.13) 



to have decay to an equilibrium state. We do not consider 
the case p = 1, because it is analogous to the previous 
one, apart fr om th e exchange of |-|-) and |— ). Therefore 
we have Eq. ( ^.11 ) and 



Pj = (^i^- 



(2.14) 



with 13 j — 2xj; by the convergence of J2j PjPj^ the com- 
plex numbers aj and Pj can be seen as the components 
of two vectors a and f3 in Z. 



Eq. (EJ) and (2.10b) give C{0)[pg] = 0; by Eqs. (t 



( ^.11 ) and ( 2.14 ) this condition reduces to [H,pg] — 0. 
Because H is s elf adj oint and defined up to a constant, 
we obtain Eq. ( ^.12aD . 
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Finally, let us choose ^ = |+). By using the rela- 
tion J2j RjRj = + P+ and by differentiating 
{N{t)) f=o two times, we obtain 

|j(iV(t))^=o + ll«|P^(iVW)/=o = 0, 

together with the initial conditions 
(7V(0))/.o = 0, 



^(iV(0)),.o 



WW 



This gives 



(^W)/=o 



WW 



t hen, condition (^.lOa ) implies f3 = Q and Eqs. ( 2.13 ) and 
( ^.14D give Eq. |2.12b| . 

The last statement of the proposition follows by direct 
computations. □ 

Now we have to find some physical restrictions on the 
possible forms of the operator S ^ U(H iS) 2). In [ pO[ , 
the case f{t) = X{t) = exp{—iu!t)9(T — t) X is consid- 
ered, where 9{x) is the usual step function and X € Z; 
for T — !■ -t-oo, X{t) represents a monochromatic coherent 
wave. Then, in 12^1 we asked 



lim lim -, — ^-^^ 

t^+ooT^+oo {N{t))l 



= 1 



VAeZ, Vw, (2.15) 



which is a form of flux conservation in the mean: if 
the possible physical processes are absorption/emission 
of single photons and direct scattering without change 
of atomic state, for large times the mean number of in- 
jected photons (A^(t))" = \\X\\t should be equal to the 
mean number of outgoing photons {N{t))\. 

The same restrictions on S arc obtained by requiring 
a balance equation on the number of photons: the mean 
number of outgoing photons up to time t plus the mean 
number of photons stored in the atom must be equal to 
the mean number of ingoing photons. 



Proposition 3 Under assumptions ( 2.1^ ), the balance 
equation 

{N{t))f + i Tr„ {a, [pit; ^ f) ~ p(0; ^, /)] } 

= {N{t)ff (2.16) 
holds yt, V^, V/ if and only if one has 

S = P+(S)S+ + P-(S)S- , S^eU{Z). (2.17) 



Proof. Any bounded operator on H (E) Z, like S, can 
always be decomposed as 



- (7_ ® F" 



(2.18) 



where 5^, are bounded linear operators on Z; the 
unitarity of S implies some simple relations among S^, 

By using Eqs. (^.5[), (^.6[), ( ^.12| ), we compute the 
t ime derivative ofTr^{CTzp(t; ^, /)}. Then, we in sert Eq. 
( ^.6b| ) into Eq. and, by using also Eq. ( |2.1^ ), we get 



{N{t))j 



+ 2Tr„{^.[p(i;^,/)-p(0;e,/)]} 



d,sTr„{[||F+/( 
+ {S+fis)\F+fis))a+ 
+ {F+fis)\S+fis))a^]p{s;C,f)} 



By the arbitrariness of t, / a nd ^, condition (2.16) is 
equivalent to F^ = and Eq. (2.17) is proved; the uni- 
tarity of follows from the unitarity of S. □ 

From now on we assume Eqs. (|Il|), (|I|), ( ^ ) to 
hold and, always for physical reasons, we take 



(2.19) 



In order to have an atom stimulated by a monochro- 
matic coherent wave we take 

f{t)^X{t) = e-''^*0{T-t)X, XeZ, w> 0.(2.20) 

The step function 6 is defined by 6{x) = 1 for a; > and 
d{x) = for a; < 0, so that X{t) represents a monochro- 
matic wave for T — s- +oo. 

Quantities like wo, ce, are phenomenological param- 
eters, or, better, they have to be computed from some 
more fundamental theory, such as some approximation 
to quantum electrodynamics. The whole model is mean- 
ingful only for uj not too "far" from wq and lvq must in- 
clude the Lamb shifts. In the final results one can admit 
a slight w-dependence in the direct scattering matrices 
5±. 



III. THE MASTER EQUATION 

In this section we study the master equation ( |2.5| ) and 
t he lo n g tirn e behav io r of the atom; the relations (^■(\), 
( ^l2| ), ( IITtI ), (^l9|) , ( |l20| ) hold. 
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The reduced statistical operator 



n = 2\{a\S-\)\ 



(3.6c) 



First of all, by setting 

p\{t) ^ lim expHa^iP + ujt)} p{t;^,X) 

T— *+oc 

xexp{~^a,{l3 + Ljt)} , (3.1) 
/3 = arg{-(5-A|a)}, 
exactly as in Proposition |l| we obtain the master equation 
d 



dt 



Px{t)^Cx [pxit)] 



(3.2) 



with the time independent Liouvillian 

C,[p]^-i[H,,p] + ^Y.{[^'P'^i\ 



At 



Rj , pRj 



5-''3(e,|a)(7_ + (ej\S+X)P+ 
+ {e,\S-X)P^, 



Hx = ^ i^o - Lu)a,^ -^\{a\S A) \ay 



(3.3a) 



(3.3b) 



(3.3c) 



The general master equation for a two-level system is 
studied in ; in the following we shall use similar tech- 
niques, apart from a different parametrization of the sta- 
tistical operator, which turns out to be more convenient 
in our case. By setting 



pxit) = 



L{t) V(t) 



v{t) l~u{t) 

< u{t) < 1 , 

u{t) > u^{t) + \v[t)\'^ 



(3.4a) 



(3.4b) 



where the conditions ( |3.4b| ) express the fact that p\ [t) is 
a statistical operator, we obtain from the master equa- 
tion 



where 



u{t) = -Gu{t) 



u{t) 
u{t) = I v{t) 




0/2 
0/2 



(3.5) 



(3.6a) 



v{t) 



6 = y ||af -i(Aa)-Im(5+A|P„S'-A)), (3.6d) 



K2 = l + ||A5A||7||af 

/\uj = oj - [ojo + Im (5'+A|Pj_S'-A)) 



(3.6e) 
(3.6f) 
(3.6g) 



5 = ^, P„ = |5)(5|, Pi = l-P„. (3.6h) 



The quantity O can be interpreted as the bare Rabi fre- 
quency. Moreover, we have 



detG=|laf (Acj)' + rV4 , (3.7) 



with 



= K^WaW^ + AK^\\afRe{S-\\Pa, {S+ + S') A) 
-4(lm(S'+A|P„S'-A))^ 
af + ||Pj_A5A|f + 2 |(5|5"A)|^ 

2Re(S'+A|Pc,S'-A)) V (S'+ + S'-)A)| 



||a||2(l + K2) + ||P^A5A|| 



(3.8) 



Let us note that ||a|| > implies det G > and > 0. 



The equilibrium state and the general solution of the master 
equation 



The equilibrium state is given by 



lim px{t) = p^ 



u{oo) v{oo) 



t — >+oo 



(3.9) 



w(oo) 1 — u{oo) 



where u(oo) and v{oo) are computed by equating to 
zero the time derivative in Eq. (3.5); then, we have 
u{oo) = G~^w, which gives 



u(oo) 



(Aa))2 +F2/4 ' 



(3.10a) 



G = 



-e"5 (a|A5A) + O 
-e-''3 (A5A|a) +0 



0/2 -0/2 \ 
b 

b J 



(3.6b) 



v{oo) 



0/2 



(Atj)2 +F2/4 



(3.10b) 



X I Yl|a|p + iAtj-|-iIm(S'+A|Pc.S'-A) 
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For the computation of the fluorescence spectrum in 
Sect ion we shall have to solve the master equation 
( |3.2[ ) also when the initial condition is not a statistical 
operator. If 



fll (712 
0'21 0-22 

is a generic 2x2 matrix, we can always write 



e"" M = (ctii + o-22)a 
where 



dsit) -di{t) 



d{t) = e-"<i(0) , d{t) 




(3.11) 



(3.12) 



lim {Yi, 

A910 \ 



a) =<5™,oi%/2[TT, (3.17a) 



lim (y,„i 



a) A ) = (5™ W2l + 1 e'^' sin Sf 



(3.17b) 



Let us recall that Yio{9, 4>) = y Pi{cos6), where the 

functions -P/(^) are the Legendre polynomials. 
Now, we set 



°° 21 + 1 ± 
= iV^^e"^' sin(5,=*=Pz(cos6') , (3.18a) 



di(0) = CTii - {an + (722)u{<x,) , 
(^2(0) = CT12 - ((Til + a'22)w(oo) , 

^3(0) = 0-21 - (0"11 + (J22)v{oo) . 



(3.13) 



Spherically symmetric atom stimulated by a collimated laser 

We end this section by particularizing our model to 
the case of a spherically symmetric atom stimulated by 
a well collimated laser. 

Let us recall that the Hilbert space Z contains the di- 
recti ons o f propagation of the electromagnetic field [see 
Eq. ( 1.10 )]. So, in order to describe a laser beam propa- 
gating along the direction 9 = 0, we have to take 

A = 7/||a|| e'*A, 77 >0, 5e[0,27r), (3.14a) 



l[o,Ae](^) 



A6iv/27r(l-cosA6') 



(3.14b) 



where l[o,Ae](^') = 1 for < 6* < A6I, l[o,Ae](^') = else- 
where; in all the physical quantities the limit A9 J, will 
be taken. Note that the power of the laser ?kj||A|p = 
hiu\\a\\^ri'^ / {AO)"^ diverges for A9 I 0, because we need a 
not vanishing atom- field interaction in the limit. 

Let us denote by Yim{9, <j)) the spherical harmonic func- 
tions; then, the spherical symmetry of the atom requires 



Y 



00 1 



1/V47r 



Im 



(3.15) 



(3.16) 



where the quantities (5/~ and are the phase shifts for 
the direct scattering in the up and down atomic states 
respectively. Let us note that we have 



A.g = g+ - 
2Auj 



s = Sg - So , 
y = z - Y sm 2s , 



■^(2Z-f l)sin2((5+-(5-) 



e = (l + v'\\P±Ag\[ 



b' 



+ 77^ (l + K2+r;2||Pj_Ag| 

1 z H sm 2s , 



(3.18b) 
(3.18c) 

(3.18d) 



(3.18e) 
(3.18f) 



/ 2 -V -V\ 
G' = 2r]&'coss b' . (3.18g) 

\ 2?7e-'''coss V J 
In order that all these quantities be finite, we require also 



^{2l + l)sin^ sf < +00, 
^(2;-f l)|sin2((5+-(5-)| < 



Then, we have 

/3 = t:-S-2So , 
bJ}^b' 



G = 



G' 



Auj — (jj ~ (wo + ifs) 
n^r,\\a\\\ 



(3.19a) 
(3.19b) 

(3.20a) 
(3.20b) 

(3.20c) 

(3.20d) 
(3.20e) 
(3.20f) 
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K2 = l + r;2||Ag|p 



u(oo) = 

v{(X)) — 



Ag{9) = i 



9 ? 

V 



47r 



sins+ (P_iA.g)(6l), 



,_i v47r 



X sin ((5;^ — (5;"") Pi(cos 
sin^ s + \\P^Ag\\^ 

CO 



|A.9|1 



1=1 



||g±||2=£(2/ + l)sin^5, 



± 



(3.20g) 
(3.20h) 

(3.20i) 
(3.20j) 



(3.20k) 
(3.201) 

(3.20m) 
(3.20n) 



1=0 



Low intensity laser. For future use, it is useful to par- 
ticularize the previous quantities to the case of a laser of 
vanishing intensity, i.e. r/ = 0: 



&' 1 - iz , 

n^o, r 



,2 

lall 



K~ = 1, C = 1 , 

^ Alu — uj — ljo ■ 



(3.21a) 
(3.21b) 



No direct scattering. The usual model of a two-level 
atom, with only absorption/emission and no direct scat- 
tering, is characterized by = 1, so that the previous 
quantities reduce to 



Alu = lo — loq . 



g± = 0, s = 0, y = z, 

= ||a 
m(oo) = 

v{oo) — 



Z2 

■n 



-I- Q'^ 



G' = 



27] 1-iz 
V 2?7 1 + iz / 



(3.22a) 
(3.22b) 
(3.22c) 

(3.22d) 
(3.22e) 

(3.22f) 



IV. DIRECT DETECTION AND TOTAL CROSS 
SECTION 

By direct detection, it is possible to measure the inten- 
sity of the light (or to count the photons) propagating in 
a small solid angle AT around some direction, which we 
take different from the direction 9 = of the incoming 
beem. The observable "number of photons in AT up to 
time t" is represented by 



7V(i; AT) 



(ei|lAT Sj) Aij{t) , 



(4.1) 



where 1at(6',(/)) = 1 for (61,0) e AT and 1at(6',0) = 
elsewhere. The fact that the direction of detection is 
different from the beam direction is expressed by 



1 



AT 



A = 0. 



(4.2) 



Then, the mean number of photons up to time t per unit 
of solid angle around (0, 0) is given by 



{n{9,cj);t)) 



1 



lAT 



■(C/(i)*(e,A)|iV(t;AT){/(i)*(e,A)) 



(4.3) 



where ^f, A(i), A are given by Eqs. (U), (|2.20D , ( |3.14| ) 
and |AT| — JJ^^ sin 9 d9d(j); the limits T +oo, 
A9 i 0, AT i {{9,(f>)} are understood. 

The (angular) differential cross section is proportional 
to the outgoing flux per unit of solid angle {n{9 , (j); t)) / 1 
divided by the incoming flux (*(^, A)|7V(i)*(^, A))/i; so 
we have 



= ^0 lim 



{n{9,4,:t)) 



t^+^ {^{t\)\N{t)^{i,\)) 



(4.4) 



where is a kinematical factor to be determined and 
with dimensions of an area. To determine Aq let us con- 
sider the cross section for direct photon scattering by 
the up or down atomic state, for which the Bohr-Peierls- 
Placzek formula (or optical theorem) gives a{9, (/>) — 
\q{9)\'^, a.j.^^ = 2^Img(0); the total cross section is 
the integral of the differential one on the whole solid an- 
gle. In our case we have to take a = and from Eq. (4.4) 

we get a{9,cl)) = |((5± - 1) A) (0, 0)| Vl|A|P and, by 
the unitarity of S^, 

^" '|(^±-l)Af 



|A|P 
2Aoy/l^A9 



Imie^'* ((5± - 1) A) (0,0) 



Then, we must have q{9) = —i^/A()A9g±{9) and, by im- 
posing the optical theorem, we get Aq = (^)^ niley ■ 
Up to now we have not taken into account the polariza- 
tion degrees of freedom. If they are taken into account 
and the cross section for not polarized light is considered, 
a 3/2 extra-factor is obtained (iQ pp. 532-533) and Eq. 
(4.4) becomes 



lim 



1 



\ Lj J 2TTri^\\a\\^ t^+oc t 



\n{9,4>;t)). (4.5) 



To compute a{9, </>) we differentiate Eq. ( [4. 3D by using 
the rules of QSC , we use Eq. ( [1.2| ) and then we apply the 
transformation (|3.lD; the final result is 
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- {n{9, 0; t)) = Tr 0)t R{9, 0) py,(t)] , (4.6) 



a-+e''r^\\a\\[g+{e)P++g-{e)P-]- 



Then, Eq. (4.5) gives 



77-^||a|pti;2 



(4.7) 



Tr{i?(0,0)ti?(0,0)p^\}. (4.8) 



Finally, by computing the trace and by using the results 
of the previous section we obtain the differential cross 
section and, by integrating it, the total one: 



Z2 + C2 



^+v'{\9+{0)?-\9-m') 



(4.9) 



47r (z2 + ^2 



■Re 



6vrc f ,1 K 

X [l + ,y2(||.9+|p-||.g_|p: 
1 



(4.10) 



Z2 + C2 



(y sin 25q + sin"' 



Let us note that the angu lar dependence in a{9, 4>) is 
entirely due to g±i 9) (3.18a ) and, so, to the presence of 
the A-term in Eq. ( |1.13D . 

By some algebraic manipulations ct^q^ can be rewrit- 
ten in a more perspicuous form: 



(^zsiuSq — cosSq)" + rj^A 



\P±9- 



||2 ^^+^ 



(4.11) 



channels are direct scattering in the up state, direct scat- 
tering in the down state and fluorescence. Some plots of 
(T^Q^ are given in Fig. 0; the independent variable 
is the "reduced" detuning z = {uj — ujQ)/\\a\\^ , the other 
parameters are given in the caption of Fig. |l|; the same 
figure contains plots of elastic and inelastic cross sections, 
which will be discussed in Sections ^ and 0. 

Whichever the line shape be, there is a strong v ariatio n 
of the c ro ss sect ion for oj around ujo+rj'^e [see Eqs. ( 3.18c ), 
( |3.18d| ), ( |3.20dD ]. The intensity dependent shift yy^e of 
the resonance frequency has received various names in 
the literature; a very suggestive one is lamp shift, a name 
suggested by A. Kastler in [Q. Note that in our two- 
level system the lamp shift is not vanishing only if the two 
states respond differently to direct scattering; moreover, 
only the contributions different from the s-wave ones do 
matter. Let us stress that also the width F of the reso- 
nance and the whole line shape are intensity dependent. 

No direct scattering. Let us also note that when the 



direc t sca ttering is negligible, i.e. when Eqs. ( 3.22 ) hold, 
Eq. ( 4.10D reduces to 



i«r/4 



G-KC^ 



w2 (Aw)2+F2/4 



(4.14) 



For a laser with negligible intensity, i.e. when rj I 0, Eq. 
(4.14) reduces to the cross section for resonant scatter- 
ing, given in |l^ pp. 530-5 33; for rj ^ 0, we have a power 
broadening (see Eq. (3.22c)) of the resonance line, which 
maintains a Lorentzian shape . 



B y com paring the general case (4.11) with the usual 
one ( 4.14 ), we see that the main differences are that in the 
general case we have lamp shift, asymmetric line shape 
and bigger power broadening. 

Lo w in tensit y lase r. For 77 = Eqs. (3.21) hold and 
Eqs. (E^gl) and ([4.1l|) reduce to 



67rc2 



9-{0) 



Gttc 



f^TOT = ll^'-L.9-| 



An {z + i) 



(2; sin Sq — cos Sq 
z2 + l 



(4.15) 



(4.16) 



A = smH+ + K^\\g+f + \\P^Ag\\ 



X l + r]' {l + \\P^Agf)sm^6o , (4.12 



V. HETERODYNE DETECTION 



B= 1 



rj'WP^AgW') (l + 7^'\\P^Ag\\ 



(4.13) 



According to the values of the various coefficients differ- 
ent line shapes appear, which are known as Fano profiles 
(Q pp. 61-63). These shapes are typical of the interfer- 
ence among various channels, when one of them has an 
amplitude with a pole near the real axis in the complex 
energy plane (see also Eq. ( 4.15| ) below); in our case the 



A. Power spectrum 

The best way to obtain the spectrum of our stimulated 
atom is by means of the balanced heterodyne detection 
scheme; the output current of the detector is represented 
by the operator |15|,p6| 



h;t) 



F{t-s)j{iy, h; ds) 



(5.1) 
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where F{t) is the detector response function, say 

J^W = fci'^exp(-It), 7>0, (5.2) 

fci ^ has the dimensions of a current, j is essentially a 
field quadrature 

jV, h; ds) = qe'"' dAh(s) + h.c. , (5.3) 

dA^,it)=J2{h\e,)dA,it), (5.4) 
j 

g is a phase factor, q Cz C, \q\ — 1, ly is the frequency 
of the local oscillator and h G Z, \\h\\ = 1; h contains 
information on the localization of the detector, say 



1 



1 



(5.5) 



where AT is again the small solid angle around {6, cj)) 
introduced in the previous section. From the canonical 
commutation relations for the fields one has 

[/(z.i,/ii;ti),/(z^2,/i2;t2)] = / dsF{ti-s) 

Jo 



xF{t2-s) (e'^'''-'->{hi\h2) 



(5.6) 



so, I{h'i, hi;ti) and /(i'2, /i2; ^2) are compatible observ- 
ables for any choice of the times either if vi — 1/2 and 
hi = /i2 either if (ft.ilft.2) = 0. Under the same conditions 
also the j's commute. 

In the following for the quantum expectation of any 
operator B we shall use the notation 

{B)l = ([/(r)*(e, A)|i3C/(T)vI/(e, A)) . (5.7) 

In the long run the output mean power is given by 



(5.8) 



^2 > has the dimensions of a resistance, it is indepen- 
dent of !/, but it can depend on the other features of the 
detection apparatus. In this section X{t) is given by Eq. 
(2.20); the limit case (3.14) will be considered in the next 



one. As a function of Piv, h) gives the power spectrum 
observed in the "channel ft,"; in the case of the choice 



(5.5) it is the spectrum observed around the direction 
{0,(f>). Proposition ^ relates P{v,h) to normal ordered 
quantum expectations of products of field operators and 
gives a sum rule which relates P{v^ h) to || A|p; let us note 
that fiwIIAlP is t he to tal power of the input monochro- 
matic state \{t) ( 2.20| ). Proposition || identifies an elastic 
and an inelastic contribution to the power and reduces 
the computation of P{v, h) to the solution of the master 
equation (3.3). For the use of QSC in the computation 
of the spectrum of a two-level atom see also Ref. |2^] . 



Proposition 4 The mean power P{i', h) can he ex- 
pressed as 



P(jy, /i) = A _|_ lim ^ 



An 2ttT 
T pt 



(5.9) 



( dAiit) dA(.)e-(*+'^)(*-^))^+c.c. , 



where k = k^k2; Eq. (5J ) holds almost everywhere in v. 
We have also 



d^= hm A(A,,(r))^ (5.10) 
I ^+00 1 



where Khh{T) = Z^ij (^l^j); moreover, for any 

c.o.n.s. {ftj} in Z, the following sum rule holds: 



E 



An 



diy = k\\X\\ 



(5.11) 



Proof. By inserting Eqs. ( |5.l| ) and (5.2) into the def- 
inition (|]^) and by changing order of integration, one 
gets 



P(u, h)^ lim — f I 



e 2 

T 



X {j{v,h;dt)j{iy,h;ds))i 



The term containing the factor exp [—7 (T ~ ^-^)] van- 
ishes for T —>■ +00 and one obtains 



P{iy,h) 



lim 



k 



T nT 



e 2 



AttT Jq Jq 
X {j{v,h;dt)j{v,h;ds))l . 



(5.12) 



By using the canonical commutation relations and nor- 
mal ordering, we have 



P{h', h) — ; — = lim 



4tt t^+oo 2ttT Jte{o,T) Jse{o,t) 
X (: j{v,h;dt)j{v,h;ds) :)l 
k 



e-*(*-^) 



lim 



e 2 



its) 



2nT Jte{o,T) 7se(o,t) 
X ({e--(*-^) dAl{t)dAhis) 

+ gV(*+^MA,,(i)dA^(s)})%c.c. 



The factor exp[ii'{t + s)], when integrated over v from Vi 
toz/2, gives rise to { exp[iz^2(i+s)] — exp[ii/i(i-|-s)] }/{i(t-|- 
s)}, which is not singular for t > and s > 0; then, the 
integ ral containing this factor vanishes for T +00 and 
Eq. (^j|) is proved. 

Now let us observe that 
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Jo Ja 

By integrating over v the second term in the r.h.s. of Eq. 
( ^.9D a Dirac delta comes out and by adding the com- 
plex conjugated term a double integral for s € (0, T) and 
t S ( 0, T) is obtained; then, by th e pr ev ious o bserv ation 
Eq. {^J§ is obtained. By Eqs. Q, < ^J^ , we 

obtain 



E 



diy^ lim -{N(T))l. 



FinaUy, by Eqs. (fj), (|^, the sum rule ^Jl\ ) 

is obtained. □ 



Proposition 5 T/ie mean power can be decomposed as 
the sum of three positive contributions 



k 

P{v, h)^ — + P^x{v, h) + Pi„ci(i^, h) , 
At: 



where 



P,i{v,h) ^ k\r{h)\ 



1 7/2 

TT {v - uSf + 7^/4 



(5.13) 



(5.14) 



-Pinol(l^, h) 



2% 



+ 00 



dt exp 



{l+\{i,-uj))t 



Tr{D(/i)t (e^^* [D{h)p^^^])] +C.C., (5.15) 



D(h) = R{h) - r{h) , (5.16a) 
r(/i)=Tr{i?(%,\}, (5.16b) 

i?(/i) = e-''3(/i|a)CT_ + {h\S+\)P+ + (/i|S'-A)P_ 

= ^(/i|e,)i?^ (5.16c) 



where R{h) is defined by Eq. (5.16c) and 



U{t) 



u{ty 



By the quantum regression theorem, which holds for a 
dynamics like U{t) [p8[, we have 



{e^''"P^{^,\)\U{t)^R{h)^U{t)U{s)^R{h)U{s)e 



X 



it A)) = Tr {i?(/i)t e^^(*-^) [R{h) e^^^ [po]] } 



where pa = exp (^(Tz/?) |^)(^| exp By recalling 

that limt^+oo e'^^*[p] = p^^ for any state p, we obtain 

lini dAlit) dAh{s)e-ii+"'y*~''>'^ 



T^+oo 27rr 







+ 00 



die 



2^ Jo 

xTr{i?(ft)te^^*[i?(/i)p,\]} 



(5.18) 



By inserting R{h) = D(h) + r{h) into Eq. ( p.l8[ ) and 
this equat io n int o Eq. ( p.9| ), we obtain the decomposi- 
tion (|5!T3|)-(|5l5|). 

The positivity of fc/(47r) and Poi(i^, ^) is apparent 
from their definitions, while to prove the positivity of 

Pincl 

{v, h) requires some transformations. 
By repeatin g in the reverse order the steps from Eq. 
to Eq. ( jtlSl ), we obtain from Eq. (|]l5|) 



Pinci{v,h) = lim — / di / dse 2 



i\t-s\ 



1 



+ lim ^ I dt dse-*l*-^l(0(s)|0(t)), 



Proof. Let us start from Eq. (5.9). We can write 

{dAiit)dA,is))l = (v|/(e,A)|C/(T)td4(t)C/(r) 
X U{T)UAh{s)U{T)^{^,X)) 

with T > t > s. By the rules of QSC (see the "output 
fields" in ||l^. Section 3), we obtain 

U{TydAhit)UiT) = U{t)''{{h\a)a- dt 

+ Y,{{h\S+e,)P+ + {h\S-e,)P.)dA,it)}Uit) . 
j 

By using this result we can write 

( dAiit) dAh{s) e^ii+'^y*-'^)^^ (5.17) 

dt / dse-(^+'(''-'^))(*~'*)(e^'"^''^'(C,A)| 
Jo 

X U{t)^R{h)^U{t)U{s)^R{h)U{s) e5 "-'^^(e. A)), 



where 



0(t) = \/ ^ c'(''-^)* U{t)^D{h)U{t)e^^'^-^{^, A) . 



By exchanging the order of integration and the names of 
the variables s and t in the second term, we get 



PineKi',/!) = hm If; I dt I dse-il*-^l 



mms)) 



which is positive because exp(— is a positive- 
definite function, i.e. the Fourier transform of a positive 
function. □ 



Notice that in the decomposition ( 5.13| ) the term 
A:/(47r), independent of v, is apparently a white noise 
contribution to the power; Pci(t/, h) is the elastic contri- 
bution, as one sees from Eq. ( 5.14] ) which gives Pe\{v, h) cx 
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d{iy — uj) for 7 i 0; fi nally , Pi^eii'-', h) is the inelastic con- 
tribution (from Eq. ( 5.15 ) one can see that no delta term 
develops for 7 J, OV 



By Eqs. (|3J), (|3.1l|)-(|3.13D, (|5.14|)-(|5.16|), we obtain 

(5.19) 



r{h) = {h\S-\) + {h\ASX)u(oo) 
+ e~''^ {h\a)v{oo) , 



/it 



27r G+'^+i{v-oj) 



d'^+c.c. , (5.20) 



2 2 



2is- « + iy 



AtT (z2 + (2) ' 



(5.26) 



(5.27) 



(/ilAS-A) 




(5.21) 



^3 



[(/i|AS'A)(l -u(oo)) 

— e~^^ {h\a)v{(X))\u{(X)) , 
[(/i|AS'A)(l -u(c5o)) 

— e~''^(/i|a)u(oo)] ^(oo) , 
e-"^{h\a) (u(c3o)- |w(oo)|2) 



(/i|AS'A)it(oo)w(cx)) . 



(5.22a) 
(5.22b) 
(5.22c) 



B. Elastic and inelastic cross sections 

Let us consider now the case of the spherically symmet- 
ric atom, sti mula te d by a well coUimated laser beam, for 
which Eqs. (|3l^)-(|3l0|) hold. We also assume that the 
dete ctor spans a small solid angle, so that h is given by 
Eq. (U) with AT i {(6',0)}, |AT| ~ sin6ld6'd0. More- 
over, we assume that the transmitted wave does not reach 
the detector, i.e. 9 > Q and so 



{h\\) =0. 



(5.23) 



From Eqs. (|1|), (|l9|)-(|^ we obtain the elastic and 
inelastic contributions to the power (per unit of solid an- 
gle) 



|AT 
1 



lATI 



P\ncl{v, h) ~ Pincl{v;0,< 



where 



P,iiiy;e,^)^k7^-'\\af\ai9)f^ 



7/(27r) 



(5.24a) 
(5.24b) 

, (5.25a) 



kr]'^\\a\\ 
+ c.c. , 



-c(0)t 



die) 

(5.25b) 



diie) 

d2{e) 

dm 



m(6') , ^ , 



(5.28a) 
(5.28b) 



+ K^y sin 2s + 2k^ cc 
+ Ag{e)K^ {k^ - iy) 



2 2 
+ iy 



(5.28c) 



mie) = Agi6) 1 - 



47r(z2 + ^2) 



(5.29) 



For the elastic and inelastic cross sections we shall have 

acl{v,0,(j)) CX Pcl{v;0,(j)), CTinol(l^; 6*, (/)) OC Pineli^; , (j)) . To 

find the co nstan t of pr op orti on ality , l et u s observe that, 
from Eqs. ( |5T0| ), ( |5^ ), (U), (|J), (U), we get 



-\-oo 



[P,i{iy;e,cf>)+Pinci{'y;e,cj,)]diy 

DO 

= hm ^{n{e,cf>;t) = &^^a{0,cf,). (5.30) 



Therefore, taking into account Eqs. (5.25), we obtain the 
expressions for the cross sections 



aei(i/;0,0) = ^|a(0)r 



7 



(zy- W)2 -f 72/4 ' 



(5.31) 



3c2 1 

crincl(i^; 0) = ^ C{0y ^ , d{0) + C.C. 

iv' G + + 1(1^ — cj) 



(5.32) 



and the relation 

f +00 



[cTciiv; 9, (j)) + (Jinc\{v] 9, (/))] dv = a{9, (j)) , (5.33) 
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where a{9, (f)) is given by Eq. (4.9). 

Finally, let us introduce the integral cross sections 

<7c\{v)^j AO sitlO j d(f)a^i{i^;e,(l)) , (5.34a) 
Jo Jo 

/ + 00 
<Jc\{iy)diy, (5.34b) 
-oo 

crinei(i^) = / desin^/ d0Crinel(i^;6', (/>) , (5.34c) 
Jo Jo 

r + oo 

Cinci = / crinci(i^)dz/; (5.34d) 



the relation ( p.33| ) becomes 

where a^^^ is given by Eqs. ( 4.12| )-(4.11) 



(5.35) 



VI. CROSS SECTIONS AND FLUORESCENCE 
SPECTRUM 

In this section we want to discuss the behavior of the 
integral cross sections and of the fluorescence spectrum. 
From Eqs. ( |5^ ), (|3^), ( |5.34a| ), ( ^.34b| ) we obtain 



7/(2^) 



(i^ - Cj)2 + 72/4 



(6.1) 



f^ci = -^^-^ ||P± [(2' + B) 5_ + Ty^AcV] I 



+ 



;t- ^ ri^K^e'* sin s — 7/ + Ik^ 
e""^o sm (5(, + — 



, (6.2) 



while from Eqs. (|5^)-(|5^), ( F^ , (|5.34cD , (|5.34dD we 

obtain 



C^2 ^2 (1 ^ ^2^ E^y^ 



67rc2 



(6.3) 



i;(2/) = (ysins + K^coss)^ + ||P_LA.g||^ (y^ + k"*) ; (6.4) 



one can check that the relation ( ^.35| ) holds true. 

Let us recall that the various quantitie s app e aring i n 
the previous formulas are given by Eqs. (pH), (|3.18a| )- 
( |3.18c| ), ( ^.18e| ), ( |3.20^ ), ( p^ and that ||apis the nat- 
ural line width, luq is the atomic resonance frequency, 
Q, = TjWaW^ is the bare Rabi frequency, rj'^e is the intensity 
dependent shift, S^, \\P_Lg±\\ , UPlA^H are parameters 
linked to the S± scattering matrices, satisfying 

P^g+\\ - WP^g^W I < |lP^A.g|l < \\P^g+\\ + \\P^g^\\ . 

We introduce also a "reduced" detuning z 



LO — UJq 



a 



(6.5) 



we have also z — 2z — 2?7^£/||q||^, s = 6q — Sq , y = 
2 

z — \ sin 2s. 

02 2 
we plot I^TOT^ f^inel, 

2 

■7^^ (Tpi as functions of the detuning z in the four cases 
rf = 10, 18, 28, 40; the other parameters are 5q = 
-0.03, 5^ = 0.13, ||P±5±ll' = 0.005, HPj^Afff = 0.02^ 
e/llajp = —0.001. Let us note the strong asymmetry in z 

of the cross sections and the fact that lim ^ cr„„„ — 

,2 



■(Tol = + sin (5q , which is about 



lim 

0.0218 with our parameters. 

Let us recall that the usual model with only the absorp- 
tion/emission process corresponds to 5^ — 0, ||^'±.g±||^ — 



A q|[ 0, £ — z = 2z; in this case, from Eqs. 

( f4.14| ), (16^), dJ), we have easily 



LO 



67rc2 ™^ 
,2 



a; 
67rc2 
2 



67rc2 



find — 



4P + 1 + 2772 ' 
42^ + 1 

(ipTTTv? 
(ipTTTv? 



(6.6a) 
(6.6b) 
(6.6c) 



Now the cross sections are symmetric in z and 



lim rr 



lim 



UJ 



■ CTcl = 0. 



2~»±oo 67rc2 ^"'^ x-^±oo 67rc2 

Then, we introduce the normalized inelastic spectrum 



67rc2 



■ O-inol(i^) 



(6.7) 



and the total one 



^TOT (''') 



67rc2 a;2 + (7/2)' 



+ Si„el(x) , (6.8) 



where we have introduced the "reduced" frequency a: and 
the "reduced" instrumental width 7 



X 



a 



7 



7 



a I 



the normalization we have chosen is 

,2 



UJ 

67rc2 



Sinoi(a;) da; = - — - CTinci • 



(6.9) 

(6.10a) 
(6.10b) 



The explicit expression of Sinci(a;) is given by the follow- 
ing proposition. 
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Proposition 6 The inelastic spectrum is given by 

.2 / 1 



Sinol(2^) 



Spinel (2;) — 



where 



1 



C = 



7r(z2 + C2)- 
+ ||P^A.9||2c"t 



le sm s 




G + 2ia; 
1 

G + 2ix 



7?||PiA5||"c' 



2 ^"t 



G' + 7 + 2ia; 
1 



d" 



c" = 



di = K m , 

+ 2k'' cos^ s + iK^ {k^ — \y) e'^ sin s , 
m' = + iy + i (z^ + - ?7^k^) e"* sin s , 

J" 2 / 2 I a2 2 2\ 



G 



(6.11) 



(6.12) 



(6.13a) 
(6.13b) 

(6.13c) 
(6.13d) 
(6.13e) 
(6.131) 
(6.13g) 



(6.14) 



G' + 7 + 2ia; 

where, by using the definitions above, 

irye'" sin s 




d + c.c. 



d = 



^2 J/ 



d = 



Then, by using the transformation 
G' + 7 = MGM-^ , 

where 

7/0 

M = I 1 



d'i 

-Tfd\ 



2^11 

3 



~7f 



/ 2 + 7 -1 ry^ \ 

2r/2e''^ cos s 6' + 7 

\-2e-'^coss 6' + 7/ 
b' is given by Eq. ^3.18^ ). 



Proof. From Eqs. ( |3^ ), ( |5^ ), ( |5^ ), ( p9| ), ( p2[ ), 
(ph, (O), (Egl), we have 



we get the Eq. ( |6.1l|) . 



□ 



By direct computations we can get the inverse of the 
matrix G + 2ij:. We can write 



1 



1 



G + 2ia; det ( G + 2ia; 



(6.15) 



where 



det ( G + 2ia; ) = (2 + 7 + 2ia;) 

+ 477^ cos s 



(k^ + 7 + 2ix) + ( z + Y sin 2s 



(k^ + 7 + 2ia;) cos s — ^-^ + ^i'^ sin s 



i?n(.T 
A2(x 
^13 (a; 



,2 ^ 2 
(k^ + 7 + 2ix)^ + ( z + Y sin2s 



+ 7 + i ( 2a; + z + Y sin 2s 



-■q 



K +7 + i(2x— z — — sin 2s 



2e coss 
2c""* coss , 



K +7 + i(2x — z — — sin 2s 



(2 + 7 + 2ia;) 



K +7 + i(2a; — z — sin 2s 



+ 2?7^e"' cos s . 



(6.16) 

(6.17a) 

(6.17b) 

(6.17c) 

(6.17d) 
(6.17e) 
(6.17f) 
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The m atrix elements D2j{x) are not needed in formula 

(in. 

One can check that the inelastic spectrum is asymmet- 
ric, but it is invariant under the transformation: x —x, 
s —s, z —z. 

The formula for the inelastic spectrum given in Propo- 
sition ^ becomes significantly simpler in the usual case 
{g± = 0) and when the intensity of the stimulating laser 
is low. 



The case gi± = 

Let us consider the usual model, when the direct scat- 
tering terms are negligible, i.e. g± — 0, which gives also 
s = 0, = 1 -I- 277^, = 1, y — z — 2z\ i n th is case the 
integral cross sections are given by Eqs. x6 and, with 
some computations, the inelastic spectrum is obtained 
from Proposition g: 



5jincl(2;) 



4?7^p(.t) 



T:q{x) (z2 + 1 + 2rff 



(6.18a) 



p{x) = (2 + 7) (1 + 7)^ + 2Tf + 



(2 + 7)^ + 2Tf + Ax"' 



(6.18b) 



+ 27 



2{2x^-Tff + {2^^f{2x^+T^^) 



If we put also 7 = 0, which means that the instrumen- 
tal width is negligible, then one can check that the flu- 
orescence spectrum Yj^^^{x), given by Eqs. (|6^), (|6^), 



( p. 181) , coincides exactly (apart from the different normal- 
ization) with the spectrum computed by MoUow ( [ p^ , 
Eq. (4.15)). Eq. ( |6.18 ) is simply the convolution of the 
inelastic part of the MoUow spectrum with a Lorcntzian 
of width 7. 

If also z = (no detuning), the eigenvalues of G can 
be c ompu ted and, by using them, the denominator in 
Eq. ( 3.1g ) can be factorized. In the case rj^ < 1/16, 

G has real eigenvalues and S]inci(a;) has a single peak 
in X = 0, while for 77^ > 1/16 two complex eigenvalues 
appear; therefore, Sinoi(a:) has a three-peaked structure 
f or ri^ sufficiently larger than 1/16. For rj very large Eq. 
( |6.18 ) gives three peaks in v ~ uj — il, v ^ oj ~ loq, 
1/ ~ cj + with height ratio 1 : ^^i^ '■ 1 and widths 

3 ll„,l|2 t „, 11^,112 I „, 3 ||„,||2 I „, -D„f [1771 „^ J^g£ 



3-1-27 
1+7 

a\\- + 7, ||q;||^ + 7; i + 7 (see Ref. [|l9|| or 
rtSf pp. 387, 423-426, 437-441 for the case 7 = 0). 



Low intensity laser 



From Eq. (6^) we see that the inelastic cross section 
vanishes in the limit of vanishing intensity of the laser; 
however, the first correction, proportional to ry^, presents 
some interesting aspects. We have immediately 



q{x) = (2 + 7) (l+7)'-|-z 



+ 4 (1 + 7) 7?2 _ 4 (4 + 3^) , 



+ Ax^ (37^ + 8J + 5 + z^ + Atj^ - Ax^) 



(6.18c) 



Now the inelastic spectrum is invariant either under the 
transformation x — > — x either under the transformation 



dire 



(4?2 + 1)^ 



Eq (z) = E{y)\^^^^ = (22;sins-|-coss)^ 
+ \\P^Agf{A^ + l) 



(6.19) 



(6.20) 



The computation of the spectrum is straightforward, but 
long; the final result is: for small i] we have 



-'incl 



(x) 



2tt 



|PLAgf(l-h7) 7(2zsins-Hcoss)^ 



1/4 



4(^2 + 1/4)-^ 



4 (a; + ?)' + (1 + 7)' 4(x-iO' + (l + 7)' 



2rf (2z sin s + cos s)^ 


'^ + {1 + 7)' /A 




7r(P-hl/4)^ 


4(x + ^' + (l+7)' 




Aix-zf + {l + ^f 



(6.21) 



In this case the inelastic spectrum is invariant either 
under the transformation x ^x either under the trans- 
formation s — !■ — s and z — > — z. 

The usual case {g± = 0) was already discussed by Mol- 
Eq. (4.30)) for 7 = and can be obtained from 
by letting vanish or from Eq. (6.21) by 
taking s = and UPlA^H = 0. In the MoUow case, for 




\Au!\ sufficiently large, the inelastic spectrum presents 
two peaks (see also Ref . p^ , pp. 106-108, 386). The 



structure given by Eq. (6.21) is similar also for s ^ 0, 
1 1 Pi A(/ 1 1 7^ 0: again two symmetric peaks appear for 
I Aw I sufficiently large. 
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Numerical computations 



In t he gen e ral ca se the total spectrum is given by Eqs. 
(|6.8|), ( 6.11 )-( 6.17 ); the analytic expression is involved, 
but plots can be easily obtained by numerical compu- 
tations. According to the values of the various param- 
eters, a well resolved triplet structure can appear, but 
also single-maximum structures can be shown. With 
the choice of parameters of Fig. Q and with an instru- 
mental width 7 = 0.6, the on resonance spectrum for 

2 

1] 



= 10, 18, 28, 40 is given in Fig. | (solid hues); the 
dashed lines give the MoUow spectrum for the same val- 
ues of and 7. The parameters in Fig. ^ have been cho- 
sen in such a way that a triplet structure appears, not too 
different from the usual one, but with a well visible asym- 
metry in the frequency x. Experiments p9|- ^ , p5| , p2[ 
confirm essentially the triplet structure; some asymme- 
try has been found, whose origin has been attributed 
to various causes. In this connection it has also been 
observed that calculations for multilevel atoms indicate 
some asymmetry Indeed, the introduction in our 

model of the interaction term containing the gauge pro- 
cess simulates the presence of other levels and the virtual 
transitions to them. 

Finally, in Fig. || we show some out of resonance spec- 
tra (detunings z = —4, —2, 3, 6) for 77^ = 28 and the 
other parameters as in Figs. |l| and ^ (solid lines); again, 
the dashed lines give the MoUow spectrum. Now, a strong 
difference from the usual case is shown, consistent with 
the strong asymmetry in z shown by the total and the 
elastic cross sections in Fig. 0. 



[1] R.L. Hudson and K.R. Parthasarathy, Commun. Math. 

Phys. 93, 301 (1984). 
[2] C.W. Gardiner and M.J. Collet, Phys. Rev. A 31, 3761 

(1985). 

[3] C.W. Gardiner, Quantum Noise (Springer, Berlin, 1991). 
[4] K.R. Parthasarathy, An Introduction to Quantum 

Stochastic Calculus (Birkhauser, Basel, 1992). 
[5] C.W. Gardiner, Phys. Rev. Lett. 56, 1917 (1986). 
[6] A. Barchielh, J. Phys. A: Math. Gen. 20, 6341 (1987). 
[7] T. Kennedy and D.F. Walls, Phys. Rev. A 37, 152 (1988). 
[8] P. Alsing, G.J. Milburn, and D.F. Walls, Phys. Rev. A 

37, 2970 (1988). 

[9] A.S. Lane, M.D. Reid, and D.F. Walls, Phys. Rev. A 38, 
788 (1988). 

[10] M.A. Marte, H. Ritsch, and D.F. Walls, Phys. Rev. A 

38, 3577 (1988). 



[11] M.J. Collet and D.F. Walls, Phys. Rev. Lett. 61, 2442 
(1988). 

[12] H.M. Wiseman and G.J. Milburn, Phys. Rev. A 49, 4110 
(1994). 

[13] A. Barchielh and G. Lupieri, J. Math. Phys. 26, 2222 
(1985). 

[14] A. Barchielh, Phys. Rev. A 34, 1642 (1986). 

[15] A. Barchielh, Quantum Opt. 2, 423 (1990). 

[16] A. Barchielh and A.M. Paganoni, Quantum Semiclass. 
Opt. 8, 133 (1996). 

[17] U. Fano, Phys. Rev. 124, 1866 (1961). 

[18] C. Cohen- Tannoudji, J. Dupont-Roc, and G. Grynberg, 
Atom-Photon Interactions: Basic Processes and Applica- 
tions (Wiley, New York, 1992). 

[19] B.R. MoUow, Phys. Rev. 188, 1969 (1969). 

[20] A. Barchielh and G. Lupieri, Photoemissive sources and 
quantum stochastic calculus, in R. Alicki, M. Bozejko, 
W.A. Majewski, Quantum Probability, Banach Center 
Publications, Vol. 43 (Polish Academy of Sciences, In- 
stitute of Mathematics, Warsawa, 1998), pp. 53-62. 

[21] H.P. Yuen, J.H. Shapiro, IEEE Trans. Inf. Theory IT-24, 
657 (1978). 

[22] A. Barchielh, On the quantum theory of direct detec- 
tion, in O. Hirota, A.S. Holevo and CM. Caves (eds.). 
Quantum communication, computing, and measurement 
(Plenum, New York, 1997) pp. 243-252. 

[23] K. Lendi, N-level systems and applications to spec- 
troscopy, in R. Alicki and K. Lendi, Quantum Dynam- 
ical Semigroups and Applications, Lect. Notes Phys. 286 
(Springer, Berlin, 1987), pp. 95-196. 

[24] A. Kastler, J. Opt. Soc. Am. 53, 902 (1963). 

[25] S. Ezekiel and F.Y. Wu, Two-level atoms in an intense 
monochromatic field: a review of recent experimental in- 
vestigations, in J.H. Eberly and P. Lambropoulos (eds.), 
Multiphoton Processes (Wiley, New York, 1978), pp. 145- 
156. 

[26] A. Barchielh, Quantum stochastic calculus, measure- 
ments continuous in time, and heterodyne detection 
in quantum optics, in H.D. Doebner, W. Scherer, F. 
Schroeck Jr. (eds.). Classical and Quantum Systems — 
Foundations and Symmetries — Proceedings of the II In- 
ternational Wigner Symposium, (World Scientific, Singa- 
pore, 1993) pp. 488-491. 

[27] H. Maassen, Rep. Math. Phys. 30, 185 (1992). 

[28] A. Frigerio, Pub. RIMS Kyoto Univ. 21, 657 (1985). 

[29] F. Schuda, C.R. Stroud Jr., M. Hercher, J. Phys. B: 
Atom. Molec. Phys. 7, L198 (1974). 

[30] W. Harting, W. Rasmussen, R. Schieder, H. Wahher, Z. 
Physik A 278, 205 (1976). 

[31] R.E. Grove, F.Y.Wu, S. Ezekiel, Phys. Rev. A 15, 227 
(1977). 

[32] J.D. Cresser, J. Hager, G. Leuchs, M. Rateike, H. 
Walther, Resonance fluorescence of atoms in strong 
monochromatic laser fields, in R. Bonifacio (ed.), Dis- 
sipative Systems in Quantum Optics, Topics in Current 
Physics Vol. 27 (Springer, Berlin, 1982), pp. 21-59. 



16 



0.06 


r| = 10 / 1 


TOT - 

inel 


0.04 


0.05 


/ ! 


el -- 




0.04 


/ ; 




0.03 


0.03 






0.02 


0.02 














0.01 


0.01 
n 




1_- 


n 



1 

- ri^= 18 


1 


1 1 1 

TOT 

\ inel 

1 el -- 




1 





-30 



-20 



-10 



10 



20 



30 



-30 -20 -10 10 20 30 



0.03 



0.02 



0.01 





FIG. 1. 



-30 -20 -10 10 20 30 -30 -20 -10 10 20 30 

3- X the integral cross sections as functions of the detuning z for 5q = —0.03, Sq = 0.13, ||PlS'±||^ = 0.005, 



IIPxAsll^ = 0.02, e/||a|| 



-0.001, and r]^ = 10, 18, 28, 40. 
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FIG. 2. Total spectrum as a function of the frequency x for z = 0, 7 = 0.6 and r/^ = 10, 18, 28, 40; solid line: 5^ = —0.03, 



: 0.13, ||Px5±lr = 0.005, IIPxApf = 0.02, £/||a|| 



-0.001; dashed line: = 0, ||Px5±|r = \\P±Ag\f = 0, £/||a|| 



0. 
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FIG. 3. Total spectrum as a function of the frequency x for 77^ = 28, 7 = 0.6 and z = —4, —2, 3, 6; solid line: S^^ = —0.03, 
So = 0.13, ||PxS±||' = 0.005, \\P±Ag\f = 0.02, e/||Q||2 = -0.001; dashed line: = 0, ||Px5±||' = HPxApf = 0, e/\\af = 0. 
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